Degree-distribution stability of scale-free networks 
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Based on the concept and techniques of first-passage probability in Markov chain theory, this 
letter provides a rigorous proof for the existence of the steady-state degree distribution of the scale- 
free network generated by the Barabasi- Albert (BA) model, and mathematically re-derives the exact 
analytic formulas of the distribution. The approach developed here is quite general, applicable to 
many other scale-free types of complex networks. 
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Introduction. The intensive study of complex networks 
is pervading all kinds of sciences today, ranging from 
physical to biological, even to social sciences. Its impact 
on modern engineering and technology is prominent and 
will be far-reaching. Typical complex networks include 
the Internet, the World Wide Web, wired and wireless 
communication networks, power grids, biological neural 
networks, social relationship networks, scientific cooper- 
ation and citation networks, and so on. Research on 
fundamental properties and dynamical features of such 
complex networks has become overwhelming. 

In the investigation of various complex networks, 
the degree distributions are always the main concerns 
because they characterize the fundamental topological 
properties of the underlying networks. 

Noticeably, for a ring-shape regular graph M of what- 
ever size, where every vertex is connected to its K 
nearest-neighboring vertices, all vertices have the same 
degree K. For the well-known Erdos-Renyi random 
graph model^ 2 ! with n vertices and to edges, the degree 
distribution of vertices is approximately Poisson with 
mean value 2m/ n. For the small- world network proposed 
by Watts and Strogatz' 1 ', the degree distribution of ver- 
tices also follows Poisson distribution approximately. 

A common feature of the above models is that the de- 
gree distribution of vertices has a characteristic size (k). 
In contrast, Barabasi and Albert^ found that for many 
real- world complex networks, e.g., the WWW, the frac- 
tion P(k) of vertices with degree k is proportional over a 
large range to a "scale- free" power-law tail: fc~ 7 , where 
7 is a constant independent of the size of the network. 
Thus, the fraction P(k) of vertices with degree k is re- 
ferred to as the degree distribution of a scale-free net- 
work. To explain this phenomenon, they proposed the 
following network-generating mechanism^ 3 ! , known as the 
BA model: 

"• • • starting with a small number (mo) of vertices, at 
every time step we add a new vertex with to (< too) 
edges that link the new vertex to to different vertices al- 
ready present in the system. To incorporate preferential 
attachment, we assume that the probability II that a new 
vertex will be connected to a vertex depends on the con- 
nectivity ki of that vertex, so that Tl(ki) = kj / y~)j hi . 
After t steps the model leads to a random network with 
t + toq vertices and mt edges." 



law form with the exponent 7 = 2.9 ± 0.1. In [4], a 
heuristic argument based on the mean-field theory led to 
an analytic solution P(k) ~ 2to 2 /c~ 3 , namely 7 = 3. To 
derive the following dynamic equation: 



dt 



mll(fcj) 
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it was assumed^ 4 ! that the probability (can be interpreted 
as a continuous rate of change of ki) for an existing ver- 
tex with degree fcj to receive a new connection from the 
new vertex is exactly equal to TO,n(fc,j), which is simul- 
taneously proportional to both the degree ki(t) of the 
existing vertex i and the number to of the new edges 
that the new vertex brings in, at time t. For notational 
convenience, this assumption will be simply referred to 
as the "TOlI-hypothesis" in this paper. 

In all the consequent works related to the BA model, 
this TOlI-hypothesis plays a fundamental role. For exam- 
ple, Krapivsky et al. ^ replaced the degree fcj (t) of vertex 
i at time t by the total number (t) of degree-A: vertices 
over the whole network at time i, thereby obtaining its 
rate equation 



dN k (t) (k - l)N k ^(t) - kN k (t) 

777, 



'km ■ 



where 5km accounts for new vertices bringing in new 
edges. In this study, the TOn-hypothesis was adopted 
in the derivations. Assuming that the steady-state de- 
gree distribution exists, using the law of large numbers 
( Nh S — > P{k) as t — > 00), they showed that the differ- 
ence equation of P(k) has an analytic solution 

P(k) 



In [3], computer simulation showed that for the BA 
model the degree distribution of the network has a power 



fc(fc + l)(fe + 2) 

for the BA model with to = 1. They also pointed out 
that only the linear preferential attachment scheme can 
lead to the scale-free structure but any nonlinear one will 
not. 

Dorogovtsev et alJ 6 ' considered fcj(t) as a random vari- 
able and defined P(k, i, t) to be the probability that ver- 
tex i has exactly k edges at time t, where vertex i is 
the vertex that was being added to the network at time 
t = i, i = 1,2, Moreover, they used the average 
of all vertex degrees as the network degree: P(k,t) = 
i ^2i—i P(k, i, t). They introduced a more general attrac- 
tion model and allowed multiple edges between vertices, 
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where each new vertex has an initial attraction degree A. 
Simultaneously, to new directed edges coming out from 
non-specified vertices are introduced with the probability 
IT, therefore k = A + q with q being the in-degree of ver- 
tices. Consequently, when every new vertex is the source 
of the to new edges like in the BA model, the attrac- 
tion model makes more sense than the BA model. They 
first arrived at the master equation of P(k, i, t) and then 
by summing all i's together they were able to derive the 
following equation: 

p (M + i) = ^-P<k-l,t)+(l-j^jP{k,t) 

To that end, by assuming the existence of P{k) [note that 
actually an additional assumption of lim^oo t[P(k,t + 
1)—P(k, t)] = is also needed], they obtained a difference 
equation for P(k). Finally, solving the equation gave an 
analytic solution 

= 2to(to+1) 
1 ; fc(fc + l)(fc + 2)' 

Here, it should be pointed out that if multiple edges are 
not allowed, then the mn-hypothesis is still needed. 

As a side note, Dorogovtsev et al.^ also considered 
the effect of accelerating growth, which is proportional 
to the power of the time variable t at each time step. 
However, this destroys the scale-free feature and degree- 
distribution stability of the network. 

Afterwards, BollobasI 8 ! made a general comment on 
the BA model: 

"From a mathematical point of view, however, the de- 
scription above, repeated in many papers, does not make 
sense. The first problem is getting started. The second 
problem is with the preferential attachment rule itself, 
and arises only for m > 2. In order to prove results 
about the BA model, one must first decide on the details 
of the model itself. It turns out to be convenient to allow 
multiple edges and loops." 

Consequently, he and his coauthors recommended a 
so-called LCD model, as follows: 

"We start with the case to = 1. Consider a fixed se- 
quence of vertices v\ , v 2 , ■ ■ ■ . We shall inductively define 
a random graph process {G\} t >o so that G\ is a directed 
graph on {vi : 1 < i < t}, as follows. Start with G\ the 
"graph" with no vertices, or with G\ the graph with one 
vertex and one loop. Given G*" 1 , form G\ by adding the 
vertex v t together with a single edge directed from v t to 
Vi, where i is chosen randomly with 

! dG i~ 1( "' ) 1 < s < t - 1 
2t_l , S 8 S I 
2t-i> s = t. 

For to > 1 wc define the process {G^} t > by 
running the process {G*} on a sequence v[, v' 2 , 
the graph G^ is formed from G™' by identifying 
the vertices v[ , v' 2 , ■ ■ ■ , v' m to form v\ , identifying 
u m+i' v 'm+2, ■ ■ ■ ."k to form w 2, and so on." 



For graph GJ^, let #„(<i) be the number of vertices 
of GJJj with in-degree equal to d, i.e., with (total) degree 
to + d, and set 

2m(m+l) 
am ' d ~ (d + m){d + m + l){d + m + 2)' 

Bollobas ct alJ 9 l rigorously proved the following result: 

lim E[#£(d)]/n = a m4 . 

Then, based on the martingale theory, they proved that 
#™j(rf)/n converges to a m .d in probability. 

It has been observed that most real-world and simu- 
lated networks follow certain rules to add or remove their 
vertices and edges, which are not entirely random. More 
importantly, at each time step, these rules are applied 
only to the previously formed network, therefore the pro- 
cess has prominent Markovian properties. Shi ct alJ 10 ' 
established a close relationship between the BA model 
and Markov chains. According to the evolution of the BA 
model, the degree ki(t) of vertex i at time t constitutes a 
nonhomogencous Markov chain as time evolves. Thus, all 
vertices together form a family of Markov chains. Con- 
sequently, based on the Markov chain theory, starting 
from an initial distribution and iteratively multiplying 
the state-transition probability matrices, the final net- 
work degree distribution can be easily obtained. Lately, 
Shi et alJ 11 ! developed an evolving network model by 
using an anti-preferential attachment mechanism, which 
can generate scale-free networks with power-law expo- 
nents varying between 1 ~ 4. There are several modified 
and generalized BA models in the literature, including 
such as the local- world BA model' 12 ! , which will not be 
listed and reviewed here. 

All in all, the BA model indeed is a breakthrough dis- 
covery with significant impact on network science today. 
Therefore, it is quite important to support the model 
with a rigorous mathematical foundation. 

It is clear from the above discussions that two key ques- 
tions need to be carefully answered for the BA model: 1) 
For the case of m > 2, can one find a scheme of adding 
new edges from the new vertex to the existing ones that 
has a probability precisely equal to mn? This is the 
key of the BA modeling. 2) Does the steady-state de- 
gree distribution of the network exist and, if so, what 
is it? This is the key to the validity of the mean-field, 
rate-equation, master-equation, and Markov-chain ap- 
proaches. The present paper will give complete answers 
to these two questions. 

Degree- distribution stability. To start, consider the 
first question. Recall that Holme and Kim' 13 ] proposed 
a scheme for new edge connection: When a new vertex 
comes into the network, the first edge connects to an ex- 
isting vertex with the preferential attachment probabil- 
ity n. After that, the rest m—1 edges randomly connect 
with probability p to the vertices in the neighborhood of 
the vertex that the first edge was connected to, or con- 
nect with probability 1 — p to those vertices that the first 
edge did not connect to. Here, consider this approach 
with p = 1 in the following scenario: When a new vertex 
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comes into the network, the first edge connects to an ex- 
isting vertex with the specified preferential attachment 
probability II, same as above. Yet, the rest m — 1 edges 
simultaneously connect to m — 1 vertices randomly cho- 
sen from inside the neighborhood of the vertex that the 
first edge was connected to. By random sampling the- 
ory this is equivalent to the above Holme-Kim scheme 
which continually connects the edges to to — 1 vertices 
randomly chosen from inside the neighborhood without 
allowing multiple edges. For this special scheme, the fol- 
lowing result can be rigorously proved. 

Proposition. For the BA model with the above special 
attachment scheme, if vertex i has degree ki(t) at time 
t, then the probability that vertex i receives a new edge 
from the new vertex at time t + 1 is exactly equal to 
mII(A;j). 

Proof. Let Pj(£+1) be the probability of vertex i receiving 
a new edge from vertex t + 1 at time t+1. Then, 



Pi(t + 1) = 



h{t) 



s~im — 2 

°fci(t)-i 



^M0~ ^ ^MO rr- J 



ieOi(t) 

E 

ieo,(t) 



y fe((t) 



m — 1 ki(t) 

T~W) =m T~W) , 



where 



C 



M*)- 



(fc,(t) - l)!/[(m - 2)!(fc,(t) - m + 1)!] _ m - 



fc,(*)V[("» -1)K*J(*) ""» + !)!] 



fci(0 



which is the probability of choosing vertex i, among the 
to — 1 vertices that were randomly chosen from inside the 
neighborhood 0/(f) of vertex I, to perform simultaneous 
connections. 

The Proposition answers the first question posted 
above and shows that the special Holmc-Kim preferen- 
tial attachment scheme is one way to implement the mlT 
hypothesis. 

In order to prove the degree-distribution stability of 
the general BA network, the BA model is specified first. 
Start with a complete graph with too vertices, which has 
a total degree Nq — mo (too — 1), and denote these ver- 
tices by —too, • • • , — 1, respectively. In all the following 
derivations, the mil- hypothesis will be assumed. The 
general BA networks will be further discussed in the last 
section below. 

Following Dorogovtsev et alJ 6 l, consider the degree 
ki(t) as a random variable, and let P(k, i, t) = P{kt(t) = 
k} be the probability of vertex i having degree k at time 
t, and moreover let the network degree distribution be 
the average over all its vertices at time t, namely, 



P(M) = 



l 



t + m 



J2 P(k,ht). 



Recall that fej(i) is a random variable for any fixed 
t and it is a nonhomogeneous Markov chain for vari- 
able i^ 10 ! . Under the mn-hypothesis, the state-transition 
probability of this Markov chain is given by 



P{ki(t + 1) = / 1 ki(t) = k} 



1 



2t+ 



2t+- 

0, 



I = 

k- 



k 
1 



(1) 



otherwise, 



where k = 1, 2, • • • , m + t — i, and i = 1, 2, ■ • • . 

The existence of the steady-state degree distribution 
for this specified BA network can be proved in three steps 
as follows. Detailed derivations are supplied in the Ap- 
pendix of the paper. 

1. Consider the first-passage probability of the Markov 
chain: 

f{k, i, t) = P{h{t) = k, h(l) ± k,l = 1, 2, • • • ,t - 1}. 

Then, the relationship between the first-passage proba- 
bility and the vertex degrees is established. 

Lemma 1. Under the mil- hypothesis, for the BA model 
with k > m, 



f(k,i,s) = P(k-l,i,s-l) 



k 



2(8 - 1) + ^ ' 

v I'm 

t-1 / 



(2) 



p(k,i,t)= n 1 



k 



s ~i-\-k—m 



2j + 



W -(3) 



2. Under the TO.n-hypothesis, using the state-transition 
probability of the Markov chain, one first finds the ex- 
pression of P(m, t), as follows: 



p(m,t) = n 1 



2i + ■ 



i + mo 
i + 1 + mo 



P(m,l) + X)- 



Z + l + m 



r n ( i 

3 = 1 V 



1 = 1 2 = 1 



. 3+ m O 

2j+Za I j + l + mo 



-!-n ■ 

: + mo 1 \ \ 



2j + 



No 



Then, one can show the existence of the limit 
limj^oo P(m, t) by using the following classical Stolz- 
Cesaro Theorem in Calculus. 

Stolz-Cesaro Theorem! 14 !. j n sequence {^}, assume 
that {y n } is a monotone increasing sequence with y n —> 
oo. If the limit lim x ™+~>-~ Xn = [ exists, where — oo < 

n^OO Vn+l-Vn ' - 

I < +oo, then lim = I. 

n — >oo y™ 

Lemma 2. Under the mn-hypothesis, for the BA model, 
the limit lim^oo P(m, t) exists and is independent of the 
initial network: 



P(m) = lim P(m, t) = - 

v ; t^oo y ' ' m + 2 



> 0. 



(4) 



3. Under the mn-hypothesis, similarly, one finds the 
expression of P(k, t) using the first-passage probability 
of the Markov chain, and then shows the existence of the 
limit limt^oo P(fc, f) by using the Stolz-Cesaro Theorem, 
if the limit limt^oo P(fc — l,t) exists. 

Lemma 3. Under the mil- hypothesis, for the BA model 
with k > m, if the limit lim t ^oo P(k — l,t) exists then 
the limit limt^oo P(k, t) also exists: 

P(k) = lim P(k, t) = P(k - 1) > . (5) 

t^oo k + 2 
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Finally, by mathematical induction, it follows from 
Lemmas 2 and 3 that the steady-state degree distribu- 
tion of the specified BA network exists. To this end, by 
solving the difference equation (5) iteratively, one arrives 
at the following conclusion. 

Theorem 1. Under the mil- hypothesis, for the BA 
model with k > to, the steady-state degree distribution 
exists, independent of the initial network, and is given by 



P(k) = 



2to(to + 1) 
k(k + l)(fc + 2) 



2TO 2 fc- 3 > 0. 



(6) 



Clearly, this degree distribution formula is consistent 
with the formula obtained by Dorogovtsev et alJ 6 l and 
Bollobas et alJ 9 l, which allow multiple edges and loops. 

Discussion. Bollobas! 8 ! once discussed the BA descrip- 
tion (the mil-hypothesis) of preferential attachment in 
detail. His result gives a range of models fitting the BA 
description with very different properties. When m > 2, 
as a new vertex comes in, it is no problem for its first 
edge to preferentially connect to an existing vertex. But 
what about the other to — 1 new edges? This question 
was not carefully addressed before. Clearly, after the 
first edge has been connected from the new vertex to an 
existing vertex, the preferential attachment probability 
n is no longer the same if later operations do not al- 
low multiple edges and loops. It is also clear that when 
to > 2, the probability of vertex i receiving a new edge 
is always greater than n. But what is it? On the other 
hand, it is also possible that the probability of vertex i 
receiving a new edge depends on other vertex degrees. 
Barabasi always emphasizes the mil-hypothesis but did 
not discuss this "how" question either. Thus, two ques- 
tions arise: 1) For the BA model, or for any other BA-like 
model, how to prove the degree-distribution stability if 
the TO,n-hypothcsis holds only approximately? 2) Is there 
a preferential attachment scheme for to > 2 such that the 
probability of vertex i receiving a new edge is indepen- 
dent of other vertex degrees? 

To answer these two questions, a new preferential at- 
tachment scheme is proposed and discussed in [15], where 



a new vertex will be simultaneously connected to to dif- 
ferent vertices and it is assumed that the preferential at- 
tachment probability n is proportional to the sum of the 
degrees - ■ ■ , k im of those vertices. They showed that 
the probability that the existing vertex i received an edge 
from the new vertex is independent of other vertex de- 
grees, namely, 



1 



too + t — m ki(t) 
to + t — 1 2mt + N m + t — 1 



where too is the number of vertices and N is the total 
degree in the initial network. Consequently, under the 
(a t ki(t) + b t )(l + o(l) fc .( t ). t )-hypothesis and some mild 
conditions, they proved the degree-distribution stability 
of Barabasi- Albert type networks. Especially, the power- 
law exponent of the network degree distribution in this 
new preferential attachment scheme is 7 = 2m + 1. 

Finally, it should be emphasized that the theory and 
scheme developed in this paper has great generality! 16 ' , 
in the sense that it can be applied to many BA-like mod- 
ified and generalized models, such as the LCD model of 
Bollobas et alJ 9 ', the attraction model of Dorogovtsev et 
alJ 6 l, the local- world BA-like model of Li and Chen! 12 !, 
and the evolving network model of Shi et alJ n l, etc. 

We summarize the results and findings in this paper 
as follows: (1) Our proving method differs from the one 
based on martingale theory, and can be applied to many 
other scale-free types of complex networks; (2) Wc do 
not need to change the BA model, e.g., to allow mul- 
tiple edges and loops; (3) We provide a special Holme- 
Kim preferential attachment scheme such that the "toH- 
hypothesis" holds. 
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To provide a rigorous proof of the degree-distribution stability of the scale- 
free network generated by the BA model, some parameters are specified as 
follows: (i) start with a complete graph with too vertices, which has the total 
degree Nq — too(too — 1), and denote these vertexes by —too, • ■ • , —1, respec- 
tively; (ii) assume that at each time step t, the probability of the new vertex 
connecting to an existing vertex i is exactly equal to mli{ki{t)). 

Here, in (ii), the preferential attachment probability is simultaneously pro- 
portional to both the degree ki (t) of the existing vertex i and the number to of 
new edges that the new vertex brings in, at time t. For notational convenience, 
this assumption will be referred to as the "TOn-hypothesis" below. 

Observe that the degree ki(t) of vertex i at time t is a random variable^. 
Let P(k, i, t) = P{ki(t) = k} denote the probability of vertex i having degree k 
at time t, and define the degree distribution of the whole network by the average 
value of probabilities of vertex degrees 

%— — mo,«7=U 

Observe also that the degree ki(t) as a process in time t is an nonhomoge- 
neous Markov chain! 10 !. Thus, for k = 1, 2, • • • , t + i — m, the state transition 
probabilities of the Markov chain, under the TOlI-hypothesis, are given by 



k 

0, otherwise. 



P{h(t + l)=l\h(t) = k}= { l = k + l (2) 



1 



1 The BA model with m = 1 



Denote the first-passage probability of the Markov chain by f(k,i,t) = 
P{ki(t) = k, ki(l) 7^ k, I = 1, 2, ■ • • , t — 1}. First, the relationship between the 
first-passage probability and the vertex degrees is established. 

Lemma 1 For k > 1, 

f(k,i,s) = P(k- l,M-l)^y^-, (3) 

Proof First, consider Eq. (3). The degree of a vertex is always nondecreas- 
ing, and increasing at most by 1 each time, according to the construction rule 
of the BA model. Thus, it follows from the Markovian properties that 

f(k,i,a) = P{ki{a) = k, h{l)^k, I = 1,2, ■■■,8-1} 

= P{h(s) = k, h{s -l) = k-l, h(l) + k, I = 1, 2, ■ ■ ■ , s - 2} 

= P{h(s) = k, ki(s-l) = k-l} 

= P{h(s - 1) = k - 1} P{h(s) = k\ki{s - 1) = k - 1} 

Second, observe that the earliest time for the degree of vertex i to reach k 
is at step k + i — 1, and the latest time to do so is at step t. After this vertex 
degree becomes k, it will not increase any more. Thus, Eq. (4) is proved. 

Lemma 2 (Stolz-Cesaro Theorem) In sequence {y-}, assume that {y n } is 
a monotone increasing sequence with y n — > oo. If the limit lim Xn + 1 ^ Xn = ; 

n^oo Vn+i Vn 

exists, where — oo < / < +oo, then lim — = I. 

n — >oo U n 

Proof This is a classical result, see [14]. 

Lemma 3 For the probability P(k, t) defined in Eq. (1), the limit limP(l,t) 

t— *oo 

exists and is independent of the initial network; moreover, 

P(l)= limP(l,t) = H> . (5) 
Proof From the construction of the BA model or Eq. (2), it follows that 

P(i,i,* + I) = (i-^)P(1,M). 

Since P(l, t + 1, t + 1) = 1, one has 

, t+i 

i— —mo ,i^0 

= ! 1 i p(i t) H 

t + 1 + m V 2t + N J v ' ; t + 1 + m ' 



2 



Then, by iteration, 



t-i 



p(i,t) = n 1 



i 



i + m 



2i + N J i + 1 + m 



P(M) + £ 



J+l+mo 



1 n(i-^ivo) 

7 = 1 



— nf 1 -^ 



2i + iV 



(i + -o)p(i,i )+ Eri(i-^- x " 



Next, let 



(i+m )p(i,i) + g n (l 2j+N ) 



and 



n-1 



4=1 ^ 

Thus, it follows that 



2i + N 



> 0. 



™ / 1 



and 



2/n+i -y n = 



3n + iV + m 
2n + AT 



2i + JV 



> 0. 



Since y„ > and ~ JM > 0, {y n } is a strictly monotone increasing nonneg- 
ative sequence, hence y n — > oo. Moreover, 



x„ + i - ar„ 2n + iV 



- (n — » oo). 



J/n+i - J/n 3n + iV + m 3 
From Lemma 2, one has 

P(l) 4 lim P(l,t) = lim ^ = lim Xn+1 - Xn = - 

t^oo n^oo y n n->oo y n+1 - y n 3 



> 0. 



This completes the proof. 

Lemma 4 For k > 1, if the limit lim P(fc — l,t) exists, then the limit 
lim P(k,t) also exists and, moreover, 



P(/c) 4 lim P(fc, t) = % \ P(k - 1) > 0. 

t— >oo k + 2 



(6) 



3 



Proof First, observe that 

P(M) - ^ ± P(k,i,t) 



i ^ . , t 1 ' 



i + m ^— ' t + m t 

i=—mo «— 1 
t 

Next, denote P(fc, t) = |^ P(k, i, t). One only needs to prove that the limit 
lim P(k,t) exists, which will imply that the limit lim P(k,t) = lim P(k,t) 

t—>oo t—>oo t—*oo 

exists. 

To show that the limit of P(k, t) exists as t — > oo, observe that P(k, i, t) = 
when i > i + 1 — fc, since in this case even if this vertex i increases its degree by 
1 each time, it cannot reach degree k. Then, it follows from Lemma 1 that 
. t+i-fe 

p(k,t) = - V P(k,i,t) 
1 i=l 

t+l-fe t t-1 / ^ 



\ e e /(*.<.-) nf 1 

i— 1 s— i+k— 1 j— s 



2.? + A^ 

., t + 1 — fe t , L _L ✓ 



' fet „SC. '2(,-l) + lV„^lV 2j + N„ 

- 7gg^- 1 ,»- 1 )^ TJ ^n(-^) 



2 — K 



t-1 7/1 1 \ ( / 7 \ "I ' 



Next, let 

(fc-1) 2 



x„ 4 P(fc-l,fc-l) 



2(fc- l) + 7Vo 



2/ + iV f_A V 2. 7 + iV 

J—K 



4 



and 



Obviously, 



y n = n J| ( 1 



i—k 



x n +i - x n = P(k - l,n) 



2i + N 



n(k - 1) 
2n + N 



> -> oo. 



n i- 



2 J + Wo 



and since 

Vn+l - Vn 



(n + l)-n(l 

(k + 2)n + A^o 



2n + 7V 



n 



i—k 



N Q + 2i 



2n + N 



n 



i — k 



2i + N a 



>0, 



one has that {y n } is a strictly monotone increasing nonnegative sequence, hence 
y„ — ► oo. Also, by assumption, 



I-n+l 



(fc - 1)n P(fc-l,n) 



1 



2/n+i ~Vn {k + 2)n + N 
Thus, it follows from Lemma 2 that 



k + 2 



lim P(k 1 t) = lim — 



lim 



/j + 2 



P(fc-l) (n->oo). 



P(fc- 1) > 0, 



therefore, lim P(k,t) exists, and Eq. (6) is thus proved. 

i — >oo 

Theorem 1 The steady-state degree distribution of the BA model with 
m = 1 exists, and is given by 



P(k) 



k(k + l){k + 2) 



~ 4k- 3 > 0. 



(7) 



Proof By mathematical induction, it follows from Lemmas 3 and 4 that 
the steady-state degree distribution of the BA model with m = 1 exists. Then, 
solving Eq. (6) iteratively, one obtains 

m = tip ( ,_ 1) = titlt^ (l _3,. 

By continuing the process till k = 3 + 1, one finally obtains 

4 



P(k) 



k(k + l)(fc + 2) 



~ 4fc~ J > 0. 



From Theorem 1, one can see that the degree distribution formula of Krapivsky 
et alJ 5 l is exact, although the mathematical proof there was not as rigorous as 
that given above. 
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2 The BA model with m > 2 



Lemma 5 Under the mil-hypothesis, the BA model when k > m satisfies 

fc-1 



f(k,i,s) 



P(k-l,i,s-l) 
t 



2(s-l) + 

p(k,i,t) = Yi n (i 

s—i+k — rn j — s \ 



No 



2.7 



No 



(8) 
(9) 



Proof First, consider Eq. (8). The degree of vertex is nondecreasing, and 
increasing at most by 1 each time, according to the construction of the BA 
model. Thus, it follows from the Markovian properties that 

f(k,i,s) = P{ki(s) = k, h(l)^k, J = l,2,--.,s-l} 

= P{ki(s) = k, h(s - 1) = k - 1, ki(l) ^k, I = 1, 2, • • • , s - 2} 

= P{ki(s) = k, h(s-l) = k-l} 

= P{ki(s - 1) = k - 1} P{ki(s) = k\h{s - 1) = k - 1} 

= P(fc-M,s-1) - fc ~ 1 N . 

2(s 1) + ^ 

Second, observe that the earliest time for the degree of vertex i to reach k 
is at step k + i — m, and the latest time to do so is at step t. After this vertex 
degree becomes k, it will not increase any more. Thus, Eq. (9) is proved. 

Lemma 6 Under the mil- hypothesis, in the BA model the limit lim P(m, t) 

t—>oo 

exists and is independent of the initial network; moreover, 



P(m) = lim P(m,t) = 

v ; t^oo y ' ' m + 2 



> 0. 



(10) 



Proof From the construction of the B A model or (2) , it follows that 



P(m,i,t+1) = 
Since P(m, t + l,t + 1) = 1, one has 



m 



P(m, i, t). 



P(m,t+1) = 



1 



t+i 



t + 1 + m 



P(m,i,t+1) 



t + mo / _ m 



f + 1 + m 
Then, iterative calculation yields 



1 



P{m,t) + 



1 



f + 1 + m 



P(m,t) 



n 



m 



2i + 



w 



i + m 
i + 1 + m 



P(m,l) + ^ 



l+l+mo 



I 



i=i n i 1 



2j+- 



j + l+rno 



G 



Next, let 



t-i i / 

m 



(l + mo)P(m,l) + ^n( 1 - 



1=1 j=i 



n-l I 



x n ±(l + m )P(m, 1) + E II ( 1 - 



i=i j=i 
and 



-l 



y n ^{n + m )Y[ > . 



i=l \ 

It follows that 

n / 

a; n+ i - z„ = I 1 

and 



2j + -^ 

• J m. 



(to + 2)n + ^ + totoq ^ 



^-^= — — nh-o^j >o- 

Consequently, 

z n+ i - z„ 2n + ^ 2 



2/n+i - J/n (to + 2)n + ^ + mm to + 2 
It follows from Lemma 2 that 



(n — ► oo). 



P(m) = lim P(m, i) = lim — = lim " +1 — = > 0. 

t->oo n->oo y„ n->oo y n+1 - y n TO + 2 

This completes the proof. 

Lemma 7 Under the toLT- hypothesis, in the BA model with k > to, if 
lim P(/c — l,t) exists, then lim P(k,t) exists and, moreover, 

P(fc) = Inn P(fc,i) = |^ip(fc-l) >0. (11) 

Proof First, observe that 

t 1 -i . t 

P(k,t) = V P(k,i,t) = V P(k,i,t)+ VP(M,t). 

v 7 t + too . ^ . Jn v ; i + too . ^ v 7 i + too t f-f 

2= — mo, 2^=0 i--mo 2=1 
t 

Denote P(k,t) ~ j ^2 P{k,i,t). One only needs to prove that the limit lim P(k,t) 



i=i 



t— >oo 



exists, which will imply that the limit lim P(k, t) = lim P(k, t) exists. 

t — >oc t — >oo 
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To show that the limit of P(k, t) exists as t — > oo, observe that P(k, i, t) = 
when i > t + m — k. Therefore, it follows from Lemma 5 that 



t+m — k 



p(k,t) = j jr p(k,i,t) 



t+m — k t t—1 I 

= !e e /(M,->n 1-577* 



i—1 s—i-\-k—m 
t-\-m — k t 



= \ E E 

z—l s— i+fc— m 
1 t s+m— fc 

= 7 E E 



fe - 1 



t-1 



V / m 3 



n ' 



fc-i 



t-i 



s— fc— m+1 i—1 
t s-1 



2(s-l) + ^ i-J 

\ I'm j=s 



7 E E^fe-Li,*-!) 



k - 1 



m 
t-1 



n > 



=fc— m+l i=l 
t 



2(«-l) + £ 



n i 



2j + ^ . 
k 

' 2j + ^ 
k \ 



2j + ^ 

J m 



= 1 t w-i,.-.)|*|ni' 



s— k— m+l 
t-1 



2( s _i) + iyo j. 



- \ n 



i—k — m+l 
t-1 



2« + ^ 

m 

i(fc-i) 



+ E ^-MJa a 

(=fe-m+l m j=fc-m+l 



P(fc- l,fc-m) 



2j + ^ i 

(fc-l)(fc-m) 

2(fe-m) + ^ 

\ -l 
fe \ 



2 ? + ^ 



Next, let 

x„ = P(fc — 1, k — m) 



(fc-l)(fc-m) 
2(fc-m) 



jVn 



E nk-W^ n 



and 



It follows that 



l—k — ra+1 

n-1 

Vn = n Y\ 

i=k—m+l 



2/ + ^ 

T m j=fc-m+l 



1 - 



2j + 



2^ + ^ 



> 0. 



x„+i - x„ = P(k - l,n) 



Tt(fc - 1) 

2n+^ . , 

m j=k-m+l 



n 



2? + ^ 

^ 777. 



and 

y«+i - = 



(n + 1) - n 1 



2n + 



iV 



n 

i—k— m+l 



2i + 



iV 
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(k + 2)n+^ 



II ( 1 2? + ) > °' 

171 + m l= k- m +l \ 11 + m / 



By assumption, 



Wn+i-i/n (fc + 2)n + ^ v ; /c + 2 1 jV ; 

It then follows from Lemma 2 that 

lim P(fc,t) = lim — = lim X " +1 - X " = ^— ^ P(k - 1) > 0. 

t _>00 „^oo y n n^oo y n+1 _ ^ /j + 2 

Thus, lim P(k,t) exists and Eq. (11) is proved. 

t— *oo 

Theorem 2 Under the mil-hypothesis, the steady-state degree distribution 
of the BA model with m > 2 exists, and is given by 

, , 2m(m +1) o , , 

p ^ = WTWTY) ~ 2 ™ fc > °- (12) 

Proof By induction, it follows from Lemmas 6 and 7 that the steady-state 
degree distribution of the BA model with m > 2 exists. Equation (11) follows 
from iteration 

/ , \ k 1 _ - _ ^ . k \ k 2/b 3 - _ 

K J k+2 K k + 2 k + 1 k y h 

till k = 3 + m. Thus, one obtains 

2m(m+l) 9 , o 

v ' k(k + l)(k + 2) 



One can see that this degree distribution formula is consistent with the 
formula obtained by Dorogovtsev et alJ 6 ' and Bollobas et alJ 9 ', which allow 
multiple edges and loops. 
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